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CSj" ■ Abstract : Recently the Sudbury Neutrino Observatory (SNO) has measured both the rates as 
well as the day-night asymmetries in the flux of neutrinos from the Sun in the charged-current (CC), 
. neutral-current (NC) and elastic scattering (ES) channels. Motivated by the small but negative day- 
q ! night asymmetry in the NC channel at SNO (1.2 standard deviations away from zero), we consider 
whether such a non-zero asymmetry can be generated in a framework where active neutrinos can 
both oscillate and decay, since the NC day-night asymmetry is zero in a frame-work that includes 
pure oscillations of active flavours. We find that the neutral current day-night asymmetry is strictly 
. zero when we include both oscillations and decay. The result holds for arbitrary number of active 
neutrino species, with some reasonable assumptions on the decay parameters. Hence, a non-zero 
asymmetry in the neutral current sector, if it survives, can only come from oscillations/decay into 
sterile flavours. We also examine the allowed region of parameter space coming from a combined 
analysis of the total event rates in the charged-current, neutral current and elastic scattering sectors 
in a two-flavour framework, and find that the neutral current rates are most sensitive to (and hence 
restrict considerably) the allowed values of the decay parameter involving the life-time of the heavier 
mass eigenstate to T2 > 10~ 4 s for eV mass neutrinos. 
PACS numbers: 14.60.Pq, 26.65. +t 
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1 Introduction 



There is now very good evidence for neutrino flavour oscillation from solar neutrino data [1, 2, 3, 4]. 
Combined analyses of these data seem to indicate [5, 6] a preference for the so-called Large Mixing 
Angle (LMA) solution in the mass-squared difference-mixing angle {5-8) parameter space (in the 
case of two-flavour mixing); on inclusion of a third flavour, these results still apply to the 1-2 
mixing angle and corresponding mass-squared difference (d~i2-uj), u = 612, since the 1-3 mixing 
angle (0 = #13) is known to be small [7].) This LMA solution is sensitive to solar matter effects 
on the propagating neutrino in a mildly energy- dependent way. The Small Mixing Angle (SMA) 
solution which is also a viable solution [5, 8] to the solar neutrino problem, on the other hand, has 
significant energy dependence due to solar matter effects. A detailed study of the various available 
data on solar neutrinos in order to isolate the allowed parameter regions in oscillation parameter 
space has been the subject of many rigourous analyses [6]. Analyses of data including decay have 
also been done for example in Refs. [9, 10, 11]. Detailed analysis of models including oscillation 
and decay, applicable to atmospheric and supernova neutrinos as well, can be found in Ref. [12]. In 
this paper, we concentrate on some constraints specifically arising from solar neutrino data (both 
from the rates and day-night asymmetries), particularly from the neutral current (NC) sector, at 
the Sudbury Neutrino Observatory (SNO) [4]. 

Information on Earth matter effects on neutrino mixing and hence flavour survival and conversion 
probabilities is available through a study of the day-night asymmetry, A ND , defined as 

aND _ q - /v • /v 
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where M % is the event rate measured during i = Night, Day. Here data from elastic scattering (ES) 
of neutrinos on electrons at Super-K [3] has been augmented by data from SNO [4] . The SNO data 
includes the day-night asymmetry in ES, in charged current (CC) v e d interaction as well as in the 
NC v x d interaction. 

The SNO and Super-K ES data are consistent with each other, within error-bars. They indicate 
a small day-night effect. The SNO CC asymmetry is positive (about 2 a away from zero) while the 
NC asymmetry is barely negative (1.2 a away from zero). 

It is well-known that oscillations among purely active neutrino flavours cannot generate a NC 
day-night asymmetry. In this paper, we examine whether it is possible to include interactions that 
generate even a small (whether positive or negative) day-night asymmetry in the NC sector. We 
find that this is impossible even with the inclusion of decay into the oscillation framework. In other 
words, if, for example, neutrino decay occurs through Majoron emission, 

vi -> V x , 

even the inclusion of both decay effects in the neutrino sector as well as of the antineutrinos that 
are produced in the decay does not yield rates that depend on the Earth-matter effects (that is, 
whether it is day or night at the time of detection). Hence the NC day-night asymmetry remains 
zero, as in the pure oscillation case. This is true for arbitrary number of (active) flavours, with 
some assumptions on the decay lifetimes. 

We then analyse the rates in the NC, CC, and ES channels in a two-flavour framework, which 
is justified since the 1-3 mixing angle is small [7]. The parameters are the two mixing parameters 
{5x2 and u), and a decay parameter a 2 that involves the life-time of the heavier mass eigenstate. It 
turns out that the contribution of the antineutrinos to the rates in any channel is small (depending 
on the life-time of u 2 ). Also, there is not much variation in the allowed oscillation parameter space 
for a large range of the decay parameter a 2 in the CC and ES channels. However, the NC channel 
is very sensitive to a 2 and thereby to the oscillation parameters. In fact, it turns out that the NC 
channel is the most sensitive measurement that constrains the decay parameter a 2 . 

In summary, the day-night NC asymmetry is strictly zero in a framework involving both neutrino 
oscillations and decay (of active neutrino species). Of course, pure oscillations cannot also generate 
a day-night NC asymmetry. The contribution of anti-neutrinos to CC and ES day-night asymmetry 
is also small. Hence a clean signature for antineutrinos and hence for decay cannot come from the 
observation of such asymmetries: the best method remains to look for such signals in the rates 
rather than in the asymmetries: this will be signalled by two neutrons in coincidence with positron 
in the CC sector at SNO. However, the decay parameter itself is rather tightly constrained from 
the NC rates data. 

In Section II we calculate the day-night asymmetries in a scenario that includes both neutrino 
oscillation and decay, for the NC, CC and ES channels. We show that the NC asymmetry is zero 
(leaving details to the Appendix), with some reasonable assumption on the decay life-times. Section 
III contains some numerical results on the allowed mixing parameter space from analysis of the CC, 
NC and ES rates at SNO [4] and the corresponding day-night asymmetries. Section IV contains a 
short discussion and conclusions. The Appendix includes some details of the day-night asymmetry 
calculations. 



2 Neutrino and antineutrino survival probabilities with de- 
cay 



We consider a framework where neutrino flavours v a are mixtures of the mass eigenstates, i^, and 
the heavier mass eigenstates v { , i ^ 1, decay with life time Tj {y\ is assumed to be stable). We have, 

i 

where the mixing matrix U a i can be expressed in terms of the mixing angles. We ignore any phase 
dependence in the mixing matrix. The time evolution of the heavier mass eigenstates includes a 
decay term apart from the phase factor, exp(— iEt): 

\uj(t)) = exp(-iEjt) exp[-t/(2 lj T J() )} \uj(0)) , 

where Tj is the life-time in the neutrino rest frame and 7j = Ej/rrij is the usual Lorentz factor. 

The (electron-type) neutrinos are produced in the core of the Sun. On their way out, they can 
interact with solar matter, and then reach the Earth. Depending on whether it is day or night 
at the detector, the neutrinos may further pass through Earth matter. The neutrino survival (or 
conversion) probabilities thus depend on the mixing matrix in the Sun, in the Earth, and the life- 
time of the heavier eigenstates. We will begin with the well-known expressions for the day-time 
probabilities. Details are included in the Appendix. 



2.1 Day rates 



The neutrino produced in the core at time to propagates as a mass eignestate z/j in the Sun until the 
time tR when it can undergo non-adiabatic transitions to another eigenstate, Vj. It then continues 
to propagate outwards, leaving the Sun at time t\ and reaching the Earth at time t 2 , where it is 
detected. The probability of a flavour a produced in the core being detected as flavour /3 on Earth 
in the day-time can be written in a factorised form [13] as, 



3 



(2) 



where the superscript S indicates the dependence on solar parameters. The braces over P indicate 
an averaging over the resonance region (or equivalently, an average over the production region [13]). 
We will always assume such an averaging and drop the braces in what follows. Note that the 
assumption of such a phase averaging restricts the neutrino mass-squared difference to be at least 
~ 10- 8 eV 2 . 

The first term in the expression above represents the probability that the produced flavour v a 
will traverse the solar matter and exit it as the mass eigenstate Vj . The second term is the (vacuum) 
probability of detection of this state as the flavour vp in the detector. The factor fP contains the 
decay information which also factors out, provided we make a single assumption about the decay 
life-times, as follows. 

From Eq. 1, we see that, apart from the phase factor that arises on propagation, the probability 
amplitude includes decay terms that involve factors of the type, 
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where (I, F) refer to any of the time intervals (t ,tn), (tn,ti), or (ti,t 2 ). Now, if the decay rate 
is very large (life-time is small compared to the time taken to traverse the Sun-Earth distance), 
most of the neutrinos would have decayed by the time they reached the Earth. This is certainly 
ruled out by observations. Hence, the relevant life-times are at least of the order of the Sun-Earth 
distance. Due to this, only the time interval (ti,t 2 ) (when the neutrino propagates from the Sun 
to the Earth) is relevant for decay, and the other decay factors involving shorter time intervals can 
be neglected. With this single assumption, the decay factor fP factorises out of the probability 
expression in Eq. 2. Note that f D = l for the case of pure oscillations as well as the lightest 
(stable) neutrino. Detailed expressions for these terms are given in the Appendix. We note here 
that fP can be written in terms of the parameter otj as jf = exp(—aij/Ej), where ctj involves the 
average Sun-Earth distance R and a matter-independent term: <x,- = Rm,j/Tj [9] where rrij are the 
(unknown) masses of the heavier eigenstates. 

Similarly, the probabilities for an electron neutrino to be detected as an antineutrino, P a ^, at the 
detector can be computed and are also dependent on the decay parameters aj. An equation similar 
to Eq. 2 holds, with the matter-dependent mixing angle, u m , being redefined with the opposite sign 
for the matter dependent term. For details, see the Appendix. 

Finally, the event rate in the day time can be obtained by multiplying the incident flux, P e p§ e 
with the corresponding cross-section for detecting up over the detector volume. 

2.2 Night rates 

It turns out that the conversion probabilities at night can also be written in a factorised form as 
the product of probabilities. While this is straightforward for the case of pure oscillations, in the 
case of decay, the factorisation again needs the assumption on the decay life-times used earlier. We 
have (see the Appendix for details), 

^ = £^?/f, (3) 
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where P^ is the same as for the day probability, but the second term includes an Earth matter 
effect that depends on the path length traversed through the Earth. Furthermore, within the 
approximations used, it turns out that = fP. Again, a similar equation (with A — > —A in 
the matter-dependent terms) holds for the case when i/g is an antineutrino. It is important that 
the factorised form of the probabilities holds in both the day and night cases, for particle and 
antiparticle production; this is a crucial ingredient in showing that the NC asymmetry is zero. 
Number conservation, or unitarity, implies that, for any j, 

E Pja = 1 j 

a 

y p E = i . 

a 

The first expression is valid in vacuum, the second in Earth. Similar expressions hold for the 
antineutrino case as well. In the neutral current channel, all flavours contribute equally. The 
expression that occurs in the NC night rates calculation below is, 

EpN \p pS fN \p pE 

ea ejJj Z^i ja i 

a j a 

= E^l/f, (4) 

j 

where we have used unitarity in the second line. This is independent of the Earth-matter effects and 
hence is the same as the dav rate exnression for the same sum that involves the vacuum probabilities. 
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where we have used fP = . This equality also holds for the sum J2 a Pea- This is a key ingredient 
in our result. 

Trivially, the following sum rule is satisfied: 

E [ P ea + Pea] = 1 • 



2.3 Day-Night Asymmetries 

We are now ready to write down expressions for the asymmetries observed at SNO. There are three 
asymmetries, corresponding to the CC, NC and ES events. We begin with the neutral current 
asymmetry. The cross-sections for the three cases are known, both for incident neutrinos and an- 
tineutrinos [14, 15]. It turns out that the cross-section for neutrinos (or antineutrinos) on deuterium 
is flavour independent in the NC channel. This is another crucial ingredient to our null result on 
the pure NC asymmetry including decay. 

The total number of events in the day-time for the NC case is given by, 



a 

= K NC $ e £ [P e D a a a (E u ) + pga^m 



where P eS involves both the energy of the produced u e and that of the decay product, V a . A similar 
expression holds for the night-time rate, with the superscript D being replaced by N everywhere. 
The overall constant factor K NC includes the detector volume (number of target atoms) as well 
as the time over which the measurement is made 1 . We now make use of the equality of the NC 
cross-sections, a a = a for all flavours, a, and write, 
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We see from Eqs. 4 and 5 that the day and night rates (when accumulated over the same period) 
are the same. Hence the NC day-night asymmetry, 

and _ n 

^-NC — u ) 

in a frame-work including both oscillation and decay of (active) neutrino flavours. This result 
does not hold if the decay occurs to a sterile flavour, since neutrinos are then "lost" and number 
conservation (the unitarity argument used above) does not hold. Finally, note that in deriving this 

1 Due to differences in lengths of day over the year, the day and night data are usually accumulated over different 
times; hence, the rates are to be normalised appropriately when calculating the asymmetry. 



result, no assumptions have been made on the details of the matter-dependence either in the Sun 
or in the Earth. 

For completeness, we list the other asymmetries in the Appendix. The CC asymmetry involves 
only the electron neutrino interactions while the ES asymmetry includes the muon neutrino as well 
(and of course the antiparticles). We now move on to a numerical analysis of the data. 

3 Numerical Calculations for rates and asymmetries 

In this section, we compute the event rates in the CC, NC, and ES channels in an oscillation frame- 
work including decay. For simplicity, we work in a 2-flavour scenario so the parameters involved are 
the mixing parameters 5 12 and u and the decay parameter ct2. In a 3-flavour frame- work, apart from 
additional mixing angles, = 9± 3 and ip = 6*23, and the corresponding mass-squared differences, 
there is another decay parameter, a 3 . However, there is very little sensitivity to the 1-3 mixing 
angle <fi [8] which is known to be small, and no dependence at all of solar neutrino rates on the 
2-3 angle ip. Furthermore, in the limit of small U e3 (and hence small 0), which is driven by the 
CHOOZ result [7], it was shown in Ref. [9] that the 0:3 dependence drops out of the expressions for 
the day-time rates so the decay is still governed only by the second-heaviest eigenstate through a 2 . 
In short, the 2-flavour analysis is sufficiently sensitive to most of the features of the model; we shall 
therefore present numerical results within such a frame-work. 

The detailed expressions for the rates in a two-flavour frame-work are listed in the Appendix. 
We shall use these formulae to explicitly compute the rates in the CC, NC, and ES channels. We use 
a chi-squared minimisation program to best-fit the free parameters and obtain the allowed regions 
in the parameter-space. We also compute, for completeness, the day-night asymmetry in these 
channels corresponding to the best-fit values of the parameters. 

3.1 Allowed parameter region from rates analysis 

As stated before, we use the cross-sections given by Ref. [15] for the neutrino-deuteron CC and NC 
cross-sections and that given by Ref. [14] for the ES cross-sections. We assume a nearly massless 
majoron so that the antineutrino produced by decay has half the energy of the initial neutrino. 

We use the fiducial volume of 0.77 ktons of heavy water and the zenith angle dependence for 128.5 
"days" of day data and 177.9 "days" of night data as given by SNO [16]. We use the normalisation 
of NC events (with kinetic energy of detected electron T e > 5 MeV) as indicated by SNO [16]. We 
also use the resolution function for the CC, ES and NC case as given by SNO. For the theoretical 
predictions, we use the standard solar model of Bahcall, Pinsonneault and Basu [17]. We use the 
Boron neutrino flux and ignore the hep contribution. Note that there are about 20% errors on the 
Boron flux of neutrinos. 

We use a simple one-slab model of the Earth with an average density of 5 gm/cc. Due to the 
latitude location of SNO, it detects neutrinos that have passed through the (outer) core of the Earth 
roughly 8.7% of the time. (SNO does not see any solar neutrinos that pass through the inner core 
of the Earth). By assuming a one-slab model that only incorporates the mantle, we have effectively 
neglected the contribution from this fraction of events. We use the 2-flavour expressions for the 
probabilities as given in the Appendix; the parameters to be determined by comparison with the 
SNO data are the mass squared difference 812, the mixing angle ui and the decay parameter a 2 . 

Before we begin the actual numerical calculation, we note that the NC event rate is amenable 
to some approximations. Since the antiparticle cross-section is small (the cross-section depends 
non-linearly on the energy and the average energy of the antiparticles is half that of the particles), 
the contribution from the antineutrino sector can be neglected. (This is also true in the CC case). 
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Table 1: Best fit solutions for fixed a 2 from a combined analysis of CC NC, and ES rates at SNO 
[4] and the corresponding x 2 P er degree of freedom. Best fit points corresponding to the so-called 
LMA and SMA regions are shown. The entire allowed parameter space is shown in Fig. 1. 
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Note that the NC events are independent of mixing in the absence of decay. Hence, when a 2 = 0, 
the second term is zero and the NC rate is independent of the mixing parameters. In fact, it is 
proportional to $cr, which is the total rate in the absence of neutrino oscillation/decay. As decay 
is turned on (so that a 2 ^ 0), the NC rate begins to depend on the mixing parameters through Pjj 
(= 1 — P^i for two-flavours). Since the NC rates are well within that predicted by the standard 
solar model, to counteract the effect of a non-zero a 2 , the parameter space in mixing gets squeezed 
in order to minimise Pjj (or to maximize P^ x ) so that the second term remains small. Since it is not 
possible to squeeze P^ beyond a point (as it is bounded by unity), larger values of a 2 get ruled out 
by the SNO result that the NC event rates are within 1 a of the standard solar model [4]. Hence 
it turns out that the NC channel is the best suited to constrain the decay parameter. Finally, such 
a squeezed solution will be a bad fit for the CC and ES rates, since here the rates are depleted 
compared to expectations. Hence a combined analysis of the CC, NC, and ES data will severely 
constrain a 2 . 

With this heuristic discussion, we proceed to calculate the rates from different channels at SNO. 
We have performed a combined chi-squared analysis of the total rates at SNO. The CC and NC 
rates are strongly anticorrelated while the CC/NC and ES rates are mildly anticorrelated; however, 
we have not included this correlation in our analysis. 

While SNO CC data favours a LOW solution (low mass squared difference, large mixing angle) 
we know that when combined with data from other experiments such as GALLEX [2], this solution 
is ruled out [8]. Hence we concentrate on solutions, including decay, around the "pure oscillation" 
LMA and SMA solutions. When a 2 = the solutions are determined by the CC and ES rates alone. 
In fact this turns out to have the lowest x 2 value when a 2 is turned on. The best-fit values to the 
parameters for different values of a 2 are shown in Table 1. It is seen that a 2 = 5 MeV is disallowed 
at 2 a. 

The 1 a and 2 a allowed parameter space from fitting the combined CC, NC and ES rates in 
(5i2,tan 2 cc>) for a 2 = 0,1,5 MeV are shown in Fig. I. There is very little difference between the 
allowed regions for a 2 = 0, 1 MeV. However, there are no solutions at 1 a for a 2 = 5 MeV. This is a 
much stronger result than the limit a 2 < 18 MeV obtained by considering only the CC and ES data 
[9]. The 2 a allowed regions are shown on the right hand side of Fig. 1. We see that the "usual" 
LMA region for a 2 = 5 MeV is disallowed by the pinching effect on the NC rates described earlier. 

This result is easy to understand if we examine the allowed regions from the individual channels. 
We find that the 1 rr allowed parameter snace in mixing' for a, eiven o-o from CC data, and ES data, 



overlap almost completely, while the allowed region from NC data is very different. The allowed 
regions coming from individual analyses of NC and CC rates are shown in Fig. 2. 

We know that the NC rate is independent of mixing parameters for a 2 = 0. Even when a 2 — 1 
MeV, the entire parameter space is allowed at 1 a level as can be seen from the left hand side of 
Fig. 2. For a 2 = 5 MeV, there is no overlap region between the allowed regions from CC/ES and 
NC data; hence there are no 1 c allowed regions for such large values of a 2 . In the LMA region, 
for instance, the allowed region from the NC channel shifts towards smaller u in an attempt to 
maximise P e ^, as discussed earlier. This results in eliminating the region in parameter space that 
is preferred by the CC rate. 

Note that we have fitted the total rates and not the energy-dependent spectrum. Inclusion of 
larger amounts of decay will distort the spectrum significantly and so it seems as if analysis of the 
spectrum may limit the decay parameters better than the rates. However, as a 2 becomes large, 
there is no allowed parameter space and hence no observable spectral distortion; a rates analysis is 
therefore unlikely to further constrain this result. This is borne out by the consistency of our results 
with the detailed spectral analysis, including decay, found in Ref. [10], where a limit a 2 < 5.7 MeV 
is obtained. 

3.2 The day-night asymmetries 

We calculate the day-night asymmetries for CC and ES channels using the best-fit values of the 
parameters for different values of a 2 as given in Table 1. We find that the CC asymmetry is roughly 
0.2 and positive for allowed a 2 in the LMA region of mixing parameters, while the ES asymmetry 
is slightly smaller and also positive. The asymmetries are much smaller in the SMA case. When a 2 
exceeds 4 MeV or so, the asymmetries are nearly zero. 

The corresponding day-night asymmetries are shown as a function of the observed kinetic energy 
of the electron in Fig. 3. There is a negligible contribution from the antiparticles, which contribute 
only at low energies, T < 5 MeV. The energy-averaged values of the SNO CC and ES asymmetries 
are shown alongside for comparison. Most of the events correspond to lower energies; hence the 
average asymmetry is closer to the value at the lower energy end. There is no significant constraint 
on the decay parameter coming from the day-night asymmetries; this of course may change with 
the availability of more data from SNO. 

4 Discussion and Conclusions 

We have shown that the neutral current day-night asymmetry is zero, in a frame- work of neutrino 
mixing that includes oscillations plus decay, with a reasonable assumption on the size of the neutrino 
life-time. We discussed the effect of decay on the neutral-current, charged- current and elastic 
scattering rates in a two-flavour frame-work. In all cases, it turns out that the contribution from 
antiparticles is small in the region of allowed decay parameter space. This is because of small cross- 
sections owing to small (roughly half) energies of antineutrinos compared to those of neutrinos. 
Hence, even if there are large decay rates, very few events will be seen. We find that the neutral- 
current rates constrain the decay parameter quite stringently. On comparing with SNO data for 
the total rates in the CC, NC, and ES channels, we get a 2 < 5 MeV, which translates to a bound 
on the life-time of the heavier mass eigenstate of r 2o > 10~ 4 s for eV mass neutrinos (equivalently, 
this is the bound on T 2o /m for an unknown neutrino mass m in eV). This is very close to the bound 
that already been obtained in a detailed analysis of solar neutrino data in Ref. [10]. This is much 
stronger than the bound a 2 < 18 MeV obtained in Ref. [9] using only the charged-current SNO 
data; in fact, it is the limiting bound that can be obtained, given the Sun-Earth distance [11]. It is 
interesting that this limit is obtained only on including the neutral-current data. 



Here we have mainly highlighted the dual role of the NC SNO data: it is most constraining of the 
decay parameter as far as total event rates are concerned; however, the NC day-night asymmetry is 
completely oblivious to the parameters of either mixing or decay. An in-depth analysis of available 
solar neutrino data, in a 2- and 3- flavour frame- work including decay, has been done in Ref. [10], 
where detailed energy correlations as well as rates from other experiments such as the Chlorine 
and Gallium experiments [1, 2], have been analysed along with the energy- dependent spectra and 
day-night asymmetries from Super-K and SNO [3, 4]. Our results are consistent with this detailed 
analysis. 

Given the allowed range of ai coming from the rates analyses, we see that the corresponding 
day-night asymmetries are not very sensitive to the mixing parameters or indeed to the decay 
parameter. Hence, the best way to put limits on, or observe neutrino decay, will be by a direct 
measurement of the rates and not the asymmetries. A distinctive signature of an "electron" detected 
in coincidence with two neutrons will signal the capture of an antineutrino. This occurs through 
the direct observation via the CC interaction in heavy water at SNO: 

I7 e d — > e + nn . (9) 

Two neutrons must be seen in coincidence with the positron. If the neutrons are not observed, then 
this process is indistinguishable from the dominant one of 

v e d — ■> e~pp , 

since water Cerenkov detectors cannot distinguish electrons from positrons. Furthermore, both the 
v e and the V e interactions have precisely the same angular dependence (which is 2 : 1 backward- 
forward asymmetric) so that angular distribution cannot be used to separate the two signals. Finally, 
even if the antineutrino energy is half that of the decaying neutrino, SNO will still be able to detect 
them through the neutron coincidence measurement. It is expected that there will soon be a 
substantial improvement in the neutron detection efficiency (with the inclusion of salt in the SNO 
detector); it may then be possible to detect these channels with good efficiency. 

Note here that while the interaction of (electron) antineutrinos on protons in water is very 
large (about 100 times the cross-section of neutrinos in water at solar neutrino energies), Super- 
Kamiokande will not see these antineutrinos since their angular distribution is roughly isotropic 
and Super-K cuts out non-forward events. 

Finally, a remark about the NC day-night asymmetry measured by SNO. After all, this was the 
motivating factor for this calculation. Currently, this asymmetry is barely negative (about 1.2a 
away from zero). If, with improved statistics, this asymmetry persists, it will be difficult to find a 
source for the effect. It is possible then that there are sterile neutrinos involved in the oscillation 
frame-work; however, the extent of conversion to steriles is limited by the total NC event rate (which 
is already consistent with the BP2000 standard solar model prediction) which admittedly has large 
theoretical errors (of about 20%). Another possible source for a non-zero NC asymmetry could be 
that the NC cross-sections for v e and v T scattering on deuteron are not the same. 

So, with the improved detection of neutrons, and the accumulation of more data, SNO will be 
in a good position to set strong bounds on neutrino decay and oscillations into sterile flavours. 

Acknowledgements: I thank M.V.N. Murthy and G. Rajasekaran for many discussions and a critical 
reading of the manuscript. 

Appendix 

We derive the neutrino conversion probabilities in a frame-work that includes both oscillation and 
decav. A flavour \u„.) that is nroduced in the core of the sun (actuallv. an electron tvne neutrino"). 



can be expressed in terms of the mass eigenstates \ui), as 
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where U ai is the usual mixing matrix that can be expressed in terms of a mixing angle uj in the 
case of two-flavours, or three angles (u, 0, tp) for the three-flavour case (we ignore any CP violating 
phases). The superscript S indicates that these are the mass eigenstates in the Sun. The effect of 
decay is to modify the phase factor that determines the time-evolution of the mass eigenstates so 
that the time evolution is given by 

\vi) t = exp[-iEit] exp[-i/(2r i )]|i/ i ) . 

The lifetime can be expressed in terms of that in the rest frame: Tj = 7« = Ei/mi. 

The neutrino, produced at a time t in the core of the Sun, propagates adiabatically upto the 
resonance point tn, where it can undergo non-adiabatic conversion to a different mass eigenstate, 
and then continues to propagate adiabatically to the edge of the Sun which it reaches at t\. Then, 
depending on whether it is night or day, the neutrino propagates upto the Earth and gets detected 
at time t%, or traverses a further distance L through the Earth and gets detected at time t 3 . 

Hence, the probability of a flavour a being detected as a flavour f3 on Earth in the day-time is 
given by 



pD 



where U is the vacuum mixing matrix, and we have averaged over the resonance region [13]. The 
decay factor fP will be defined below. Here represents the non-adiabatic transition probability 
amplitude from i — > j so that |M^| 2 is the usual Landau Zener transition probability [18] in the 
Sun. 

As has been pointed out earlier [13], it is possible to express the probability as a factor that is 
solar- matter dependent, times a solar- matter independent factor: 



pD 
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where 



p s . = V \U S \ 2 \M 
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ji\ 



and 



As explained in the main part of the text, the decay factor fP in the expression for the probability 
amplitude involves a product of terms: 



~ exp 
= exp 



F dt_ 

F rrii dt 



-I 

~ Ii 2- , E 



with one such term coming from every propagation from tj to tp of a neutrino z/j. In general, 
we can assume the modification of the energy eigenvalue due to the matter-dependent effects is 
small, so that the integral above reduces to ~ {tp — ti)[m i /(2To i E i )] . If the time interval of 
propagation is short, then the numerator is small compared to the life-time unless the latter is itself 
very small. In the latter case, we would expect most of the produced neutrinos to have decayed 
before reach in or the Earth, but this is clearlv not true. We therefore assume that the life-time of the 



heavier mass eigenstate is at least the Sun-Earth distance (divided by the speed of light). Then the 
decay exponentials are all small and can be dropped, except the terms containing the time interval 
(tp — tj) = (t2 — ti) which is the Sun-Earth separation. The relevant integral contains the vacuum 
energy eigenvalue and hence is anyway matter-independent. The decay term that appears in the 
expressions for the probability is then the square of the above, with no cross terms, for the energy 
eigenvalue Ej corresponding to Vj that propagates in vacuum from the Sun to the Earth, 



fj = expi-aj/Ej) , 

where we have used aij = Rmj/T Q . [9] where R is the Sun-Earth distance. 

A similar formula holds at night-time, except that the vacuum probability Pjp is replaced by an 
earth-matter dependent one: 

pN _ \^ pS pE rN f \ 
r aB — / , r ai r i8J j j \ z ) 

j 

with the solar probability the same as before, and with 



pE 



E UfiU^M*Mfi exp \-i [ t3 (E* - E§)dt 
kk' i Jt 2 



Here Mj^ = (i^^Wj) is the transition amplitude from j — > k, where is the mass eigenstate just 
inside the Earth. We have assumed only one non-adiabatic transition occurs, when the neutrino 
encounters the discontinuity in density at the Earth's surface. However, it is straightforward to 
generalise this to an arbitrary number of discontinuous jumps in Earth's density, using the so-called 
slab model of the Earth. Details may be found in Ref. [13]. The decay factor is the same as in 
the day-time case since the matter effects (whether in Earth or in the Sun) due to decay are small 
due to the small distance of propagation in these. So the vacuum contribution is the only one that 
remains here as well. So we have 

ff //' • (3) 

Finally, the event rate can be obtained by multiplying the incident flux, P e ^$ e with the correspond- 
ing cross-section for detecting v@ over the detector volume. 

Similar formulae hold for antineutrinos produced during decay, with all matter-dependent terms 
being redefined with A — > —A, where A is the matter dependent term in the Sun or Earth. It is 
important that the factorised form of the probabilities holds in both the day and night cases, for 
particle and antiparticle production; this is a crucial ingredient in showing that the NC asymmetry 
is zero. 

Number conservation, or unitarity, implies that, j, 

a 

yp E = i . 

a 

The first expression is valid in vacuum, the second in Earth. Note that the second equation is 
valid even when Eq. 2 for the night-time probability is replaced by a more complicated one using a 
many-slab model of the Earth. 

Similar expressions hold for the antineutrino case as well. In the neutral current channel, all 
flavours contribute equally. The expression that occurs in the NC night rates calculation is, 



EpN _ pS fN pE 

ea / j ejJ j / j in 1 

a j a 

= E^/f, (4) 



where we have used unitarity in the second line. This is independent of the Earth-matter effects and 
hence is the same as the day rate expression for the same sum that involves the vacuum probabilities, 
P • 

a j a 

= EPeW, (5) 
3 

where we have used = f^. This equality also holds for the sum J2 a Pea- From this it is clear 
that the NC rates are the same during the day or night for arbitrary numbers of active neutrino 
species; hence the NC day-night asymmetry is zero as explained in the main part of the text. If the 
decay is into a sterile flavour, of course, this is no longer true. 

We now specialise to the two-flavour case where we list the detailed expressions for the proba- 
bilities needed for the computation of the rates and hence the asymmetries in the different channels 
at SNO. Then with the choice m 2 > mi, v\ is stable so that r 2 is the only relevant decay parameter. 
We have, 

f 2 D = f 2 N = exp[-a 2 /P 2 ] . 

We now list the various probabilities. Since the Earth's density is much smaller than that in the 
solar core, there is only a small modification of the energy due to Earth-matter effects. Then we 
can define 

X = (t 3 -t 2 )6g/(2E) , (6) 
where 5f 2 is the mass-squared difference in (Earth) matter: 

Sf 2 = 5 12 cos(2(u; — oo e )) — A E cos 2uj e , 

where the Earth matter dependent term is, Ae = 2y/2G E N^ E with Nf the electron number density 
in the Earth and (t 3 — t 2 ) = 2R E cos0 z and the matter mixing angle is u e . Here 2R E is the diameter 
of the Earth (= 5.1i?£(m)MeV/eV 2 ) and 6 Z is the zenith angle of propagation of the neutrino with 
respect to the instantaneous direction of the Sun. 

The Day-time probabilities: In the two-flavour case, P e2 = 1 — P el ; furthermore, we have, 

PZ = c 2 pS 1 + s 2 Pgexp(-a 2 /E 2 ) , 

K = s 2 P e s l + c 2 Pgexp(-a 2 /E 2 ) , 

P£ = c 2 P/ 2 [l-exp(-a 2 /£ 2 )] , 

P% = s 2 P e s 2 [l-exp(-a 2 /E 2 )}. (7) 

Here c and s refer to cos uo and sincu respectively. Defining the Landau-Zener probability P E z [18] for 
transition between the v x and v 2 eigenstates (note that the third mass eigenstate, even if included, 
does not participate in non-adiabatic transitions in the Sun, provided the relevant mass squared 
difference is that which solves the so-called atmospheric neutrino problem), we have 

P? 1 = s 2 m P LZ + c 2 m (l-P LZ ) , (8) 

where s m and c m refer to sin uj m and cos u m in the solar matter, given by 

5 12 cos2u-A 

cos 2u m = . (9) 

'(5f 2 sin 2 2u + (5 12 cos 2u - A) 2 ) 

The above formula holds for antineutrinos with A — > —A and for neutrinos and antineutrinos in 
Earth matter as well, with A — > At? so that uj m (+A) — > d„(-\-A tA. 



Note that 
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Pg + P e s 2 exp(-a 2 /E 2 ) , 
exp(-oi 2 /P 2 )) , 



f 2 ( 



(10) 
= 



so that the total conversion probability into neutrinos and antineutrinos is conserved. When a 2 
we recover the usual formulae for pure oscillation case. 

Night-time probabilities: Here we need to compute, in addition, the Earth-matter effects. We have 



P E - 1 

r \e — 1 
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pE _ 



2e 



s cos 2u; e + 2csc e s e cos2a; e 
+2c e s e (c e s e cos 2u> — cs cos 2u; e ) cos X , 
1 — 2c 2 s 2 — c 2 cos 2 2cj e — 2csc e s e cos 2u; e 
— 2c e s e (c e s e cos 2u — cs cos 2u; e ) cos X , 
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(11) 



where X has been defined above and oj e is the Earth- matter dependent mixing angle. Notice that, 
along with the usual number conservation, Pf e + Pf = 1, we have also Pf e + P 2e = 1 . Substituting 
in the transition probabilities, P^, this gives us, 



-,7V 



P 



A? 



P^P* + P*) + Pg(P 2 E e + Pg) exp(-a 2 /E 2 ) 
Pf 1 + (l-Pj>)exp(-a 2 /P 2 ) . 



(12) 



Notice that the sum is independent of Earth matter effects and is exactly equal to the "Day" 
formula; this was already shown in the general case of arbitrary number of flavours. 

For the detection of antiparticles at night, we have similar expressions to those for the neutrinos, 
with uj e (A E ) — >• uj e (-A E ): 



P N -+P 

ee 1 e/i 



N 



P^Pg + P(- 

(l-PfJ(l-exp(-a 2 /P 2 )) . 



Again, the total probability is conserved and is independent of the Earth-matter effects. 
The event rates: With these probabilities, we now define the rates for CC, NC and ES events at 
SNO. We begin with the neutral current channel. The event rate at SNO in the NC channel is 
given by 

j\f NC = A/" 6 + A/ 7 * + + A/ 77 . 



Hence, 



Af NC 

J *Day 



■ /v Night 
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(13) 
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exp(-a 2 /P 2 )) (cr - cr) 
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(14) 



These expressions leads us to the general result that we had shown earlier: that the NC day-night 
asymmetry is zero. Here a (a) is the flavour-independent neutral current cross-section of neutrinos 
(antineutrinos) on deuteron. For calculating the cross-section and hence the rates, we make the 
simnle assumption that in the decav. the antinarticle carries awav about half the enerev nf the narent 



neutrino. Hence, although the NC cross-sections for neutrinos and antineutrinos are roughly the 
same at the same energy, the contribution from antineutrinos is very small since the cross-section 
rises exponentially with energy Indeed, for most values of the bulk of the contribution to the 
rates is still from neutrinos. 

Similar expressions can be written for the CC and ES rates. In the former, only the v e (or V e ) 
contribute. We have, 



• /v Day 

Af cc 

Jy Night 



'cc 

pV 4- P N rr B 



(15) 



where the probabilities involved have been listed earlier and the cross-sections are taken from 
Ref. [15]. There is no simplification of the expressions although the anti-neutrino contribution is 
small here as well. 

The ES rates involves all particles and antiparticles, but with flavour- dependent cross-sections 
as given in Ref. [14]. We have, 



KfES 
JV Day 
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P D rr e 4- P D rr^ 4- P D n e 4- P D rr^ 

' P N n e 4- P N n-V 4- P N ' n* 4- P N ' rF 
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(16) 



The cross-section for elastic scattering of the v e is about 6 times larger than than of (and even 
larger than the antineutrino cross-sections), although the values are energy dependent. 
The day-night asymmetries: The CC and ES asymmetries are not zero; the former depends only on 
P ee and Pge. These are different for the day and night-time. Hence, the CC asymmetry is non-zero 
by an amount proportional to 



(17) 



where i — e, e. It is straightforward to similarly write an expression for the ES asymmetry, which 
has both charged- current and neutral-current contributions: 



^t - XE!y = $E [(Pei ~ ^)^5 



(18) 



where the summation includes contributions from electrons and muons (and their antiparticles as 
well). 
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Figure 1: 1 a and 2 o allowed regions in (<5i2-tan 2 uS) space from a combined analysis of CC, NC 
and ES rates at SNO [4] for fixed values of the decay parameter a 2 = 0,1,5 MeV (dotted, dashed 
and solid lines respectively). Notice that there are no 1 a allowed regions for a 2 = 5 MeV. The 
best-fit values (that correspond to a 2 = 0) of the parameters in the LMA and SMA region as given 
in Table 1 are shown as solid dots in the figure. 
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Figure 2: 1 a allowed regions in (5 12 -tan 2 cj) space from individual CC and NC rates at SNO 
(calculated with respect to the combined minimum as listed in Table 1) for different values of the 
decay parameter a 2 . The NC allowed regions are shaded for contrast. Notice that there is no 
overlap between CC and NC allowed regions for a 2 = 5 MeV. The allowed regions from ES rates 
almost completely overlap the CC ones and hence are not indicated separately. 




Figure 3: The day-night asymmetries for CC and ES (solid and dashed lines) rates for the best fit 
parameters for different values of a 2 as a function of the kinetic energy of the observed electron. The 
lines correspond to a 2 — 0, 1, 5 MeV with decreasing magnitude of asymmetry. The energy-averaged 
asymmetries measured by SNO are shown on the right hand side of the figure, for comparison. 



